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Motivated by the observation of localized traveling-wave 
states ('pulses') in convection in binary liquid mixtures, the 
interaction of fronts is investigated in a real Ginzburg-Landau 
equation which is coupled to a mean field. In that system 
the Ginzburg-Landau equation describes the traveling-wave 
amplitude and the mean field corrsponds to a concentration 
mode which arises due to the slowness of mass diffusion. For 
single fronts the mean field can lead to a hysteretic transition 
between slow and fast fronts. Its contribution to the interac- 
tion between fronts can be attractive as well as repulsive and 
depends strongly on their direction of propagation. Thus, 
the concentration mode leads to a new localization mecha- 
nism, which does not require any dispersion in contrast to 
that operating in the nonlinear Schrodinger equation. Based 
on this mechanism alone, pairs of fronts in binary-mixture 
convection are expected to form stable pulses if they travel 
backward, i.e. opposite to the phase velocity. For positive 
velocities the interaction becomes attractive and destabilizes 
the pulses. These results are in qualitative agreement with 
recent experiments. Since the new mechanism is very robust 
it is expected to be relevant in other systems as well in which 
a wave is coupled to a mean field. 
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I. INTRODUCTION 

Among the most striking recent findings in the investi- 
gation of pattern-forming dissipative systems have been 
the observations of stable localized structures surrounded 
either by the structureless basic state or by other struc- 
tures. Steady structures of this kind have been found 
in experiments on Rayleigh-Benard convection in narrow 
channels where patches of rolls with longwavelength co- 
exist with patches of short wavelength ]l]~|| . In Taylor 
vortex flow large, turbulent axisymmetric vortices have 
been found to coexist stably with small, laminar vor- 
tices [0. Another class of localized structures is given 
by the localized drift waves observed in a variety of sys- 
tems undergoing a parity-breaking bifurcation. They 
were first found in directional solidification ||. Subse- 
quently they were also observed in other interface in- 
stabilities like the printer instability f|] and in cellular 



flames 0. Based on theoretical predictions ||, the same 
phenomena were then found in Taylor vortex flow be- 
tween counter-rotating cylinders Q . Theoretically, these 
localized drift waves have been successfully described 
with phase-amplitude equations adequate for a parity- 
breaking bifurcation ||,|l(J-|l6| . 

The present paper is motivated by localized waves 
which have been observed in the convection of bi- 
nary mixtures. They consist of traveling- wave packets 
('pulses') which drift through the conductive, structure- 
less state. They have been investigated in great detail 
and presumably constitute the best studied structures of 
this kind Jl7|-pT|. Depending on parameters, two seem- 
ingly different classes of pulses have been found. For 
weakly negative separation ratio ip of the mixture the 
pulses comprise only a few wavelengths and their size is 
relatively independent of the parameters (pulses of 'fixed 
width'). For more negative values of the separation ratio 
long packets have been observed. In early experiments 
their width varied from run to run in a seemingly ran- 
dom manner (pulses of 'arbitary width'). In those ex- 
periments both types of pulses were stationary, with the 
convection rolls traveling through the pulse. In numerical 
simulations of the Navier-Stokes equations (2^j2^| and in 
subsequent improved experiments it was found that the 
short pulses [^0| and then also the long pulses in fact 
drift as had been expected on general symmetry grounds. 
The drift velocity was, however, found to be extremely 
small. 

Analytically, these pulses have been described within 
the framework of a complex Ginzburg-Landau equation. 
Two approaches have been taken. In the strongly disper- 
sive limit the Ginzburg-Landau equation can be viewed 
as a dissipatively perturbed nonlinear Schrodinger equa- 
tion and it has been shown that the pulses can arise from 
the solitons of that equation p5|~p8[ . For weak dispersion 
localized waves have been described by fronts which can 
form stably bound pairs due to the repulsive interaction 
arising from the dispersion p9| , pfj| . In either case, the 
localized waves should travel with a velocity which is es- 
sentially given by the linear group velocity of the waves. 
The experimentally and numerically observed velocity is, 
however, by a factor of 20 to 30 smaller p0j , p2|| . In addi- 
tion, the numerical simulations revealed a striking behav- 
ior of the concentration field, which is advected by the 
traveling rolls. It was suggested that it could be respon- 
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sible for the slow-down of the pulses \22\. This finding 
motivated a detailed analysis of the derivation of the con- 
ventional Ginzburg-Landau equations. It was shown that 
these equations break down already quite close to thresh- 
old; due to the slow mass diffusion in liquids, which is 
characterized by their small Lewis number L = O(10 -2 ), 
certain modes of the concentration decay very slowly and 
have to be treated as independent dynamical variables in 
addition to the convective amplitude A pl|-|33|| . For free- 
slip-permeable boundary conditions extended Ginzburg- 
Landau equations for A and a large-scale concentration 
mode C have been derived. In contrast to the conven- 
tional Ginzburg-Landau equation, they do not become 
singular in the limit £ -» and constitute a minimal 
model to describe this important physical effect. 

Numerical simulations of the extended Ginzburg- 
Landau equations in the strongly dispersive case showed 
that the concentration field can indeed slow down the 
pulses considerably |3l|-|33fl ■ Analytical investigations us- 
ing soliton perturbation theory confirm this result |53],[34j . 
Numerically, it has also been found that (long) pulses can 
arise stably even without dispersion when the concentra- 
tion mode is included pljjffi] . This shows, that in this 
system a second mechanism, which is unrelated to that 
responsible for the solitons of the nonlinear Schrodingcr 
equation, may be able to induce the localization of waves. 
In the present communication we discuss long pulses con- 
sidering them as a bound pair of fronts. We isolate their 
localization mechanism analytically and show that the 
interaction between fronts can become repulsive if the 
pulse travels backward, i.e. opposite to the velocity of 
the waves. This mechanism is quite robust and may be 
expected to be relevant in other systems as well. Equa- 
tions similar to those discussed here have been derived, 
for instance, for two-fluid channel flow j3f|. 

Recently, the interaction of fronts in a system of 
reaction-diffusion type with two fields has been inves- 
tigated (3^). In that work stably bound pairs of fronts 
have been found as well. Their localization mechanism is, 
however, quite different from that discussed here; while 
in those equations the coupling of the two fields arises 
through reaction terms, it involves gradient terms in the 
extended Ginzburg-Landau equations discussed here. 

The organization of the paper is as follows. In sec. II 
the extended Ginburg-Landau equations are discussed. 
The equations describing the interaction of fronts are de- 
rived in sec. III. They are analyzed in detail in sec. IV. 
In sec.V the effect of terms neglected in the derivation 
of the equations describing the interaction of the fronts 
is discussed using numerical simulations of the extended 
Ginzburg-Landau equations. In the concluding sec. VI we 
comment on the relevance of the concentration mode in 
combination with dispersion for understanding very re- 
cent experiments |2l|. 



II. 



THE EXTENDED GINZBURG-LANDAU 
EQUATIONS 



As shown in 31 33] the derivation of the conven- 



tional Ginzburg-Landau equations from the Navier- 
Stokes equations becomes invalid in the limit of vanish- 
ing mass diffusion of the second component. Since the 
Lewis number C, which is given by the ratio of mass dif- 
fusion to heat diffusion, is O(10 -2 ) in liquids this implies 
that the conventional Ginzburg-Landau equation yields 
a poor description of the experimental system already 
for quite small convection amplitudes. Technically, in 
the limit C — > a certain concentration mode C cannot 
be eliminated adiabatically in favor of the traveling- wave 
amplitude A and has to be kept as an independent dy- 
namical variable. For free-slip-permeable boundary con- 
ditions this leads to the following equations as a minimal 
model to capture the dynamics of the concentration mode 
& 

d- t A + SdxA = dd?A + aA + £\A\ 2 A + p|i| 4 i + fCA, (1) 
d~ t C = d c d\C + a c C + hd$\A\ 2 + gC\A\ 2 , (2) 

where the amplitude A denotes the traveling-wave am- 
plitude and is identical to that appearing in the conven- 
tional Ginzburg-Landau equation. The streamfunction 
tj) and the deviations from the linear temperature profile 
as well as that from the concentration profile are there- 
fore given by (tp,8,rj) = e Ae lqx+lut sin7rz(V>o, #0, 'To)- 
The new concentration mode C characterizes a long- 
wavelength mode of the concentration field which de- 
pends only on the vertical coordinate z, ("0, ff) = 
eC sin 2-7rz(0, 0, 1). From a more general perspective, 
eqs.(0^) describe the interaction of a finite- wavenumber 
traveling wave with a zero-wavenumber steady mode. 
The interaction of two steady modes with very different 
wave numbers has been treated in [ [37f . 

To be precise, eqs.([j],|]) do not contain all terms up to 
fifth order. The complete set up to third order includ- 
ing the values of the coefficients^ for free-slip-permeable 
boundary conditions have been presented in | p2| . Here 
we want to focus on the advection of C by the waves 
and the resulting feed-back on the traveling-wave ampli- 
tude through the dependence of its growth rate on the 
local (vertical) gradient in the concentration field. These 
important effects are described in eqs.([j],|]) by the terms 
/i9j|A| 2 and fCA, respectively. In addition, the damp- 
ing term a c C (a c < 0) arises from vertical diffusion of the 
concentration. Through the term C|A| 2 this damping is 
also affected by convection. Note that eqs.(|l],|]) are only 
adequate if the convective amplitude is small enough to 



l In jB2| it had been overlooked that the coefficient /14 has to 
vanish in general due to the fact that the concentration mode 
C is real IBSf. 
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satisfy a c + g\A\ 2 < 0. If this term were positive, plain 
waves \A \ 2 = A 2 , would be unstable to an evergrowing 
concentration mode. The leading-order contribution to 
the diffusive term d 2 C does not arise from diffusion of 
the concentration; instead large-scale variations in the 
concentration field affect the local bouyancy of the fluid 
and generate vorticity which in turn advects the basic 
(linear) concentration profile. 

To focus on the new localization mechanism, which 
was found numerically not to rely on the dispersive ef- 
fects embodied in the imaginary parts of the coefficients 
p2| , we consider here the special case in which all the 
coefficients and the amplitude A are real. In addition, to 
study the dynamics and interaction of fronts we restrict 
ourselves to the regime below threshold, a < 0. Eqs.(0J|) 
can then be simplified by introducing scaled variables via 



~\ i/4 r 

?=) A, C = -lc, 
a I a 




a = — , h 
a 



c d c 
lap d 



hf 



-pa 2 d 



'ap 



This leads to 



d t A + sd x A = rj 2 d*A - A - 
d t C = 5dlC - aC 



cA 3 -A b + CA, 
- hd x A 2 + gCA 2 . 



(3) 
(4) 
(5) 
(6) 



(7) 

(8) 



Note that in the scaling (|[^) the reduced Rayleigh num- 
ber a has been scaled to -1. The parameter which allows 
a scanning between the Hopf bifurcation at a — and the 
saddle-node bifurcation at a = c 2 /Ap is now given by c. 
Its values range therefore from c = 2 at the saddle node to 
c — ► oo at the Hopf bifurcation. For free-slip-permeable 
boundary conditions the coefficients a, <5, g and, in par- 
ticular, h are positive J32| . In ([I|) we have introduced the 
parameter rj. It characterizes the strength of the diffusion 
of the traveling-wave amplitude A and therefore governs 
the width of the fronts in A. In this paper we will con- 
sider the limit of narrow fronts, i.e. i) C 1. This allows 
the derivation of equations describing the dynamics and 
interaction of fronts which focusses on the contribution 
from the concentration mode. 



III. DERIVATION OF THE FRONT EQUATIONS 

The goal of this paper is an analytical description of 
the contribution of the concentration mode to the inter- 
action of fronts. In the standard approach to capture the 
interaction of fronts analytically one assumes the distance 
between the fronts to be large. The interaction, which is 



essentially due to the overlap of the fronts in the convec- 
tive amplitude, is then weak and the resulting dynamics 
slow. In the present system the fronts interact addition- 
ally via the concentration mode. In certain parameter 
regimes the decay length of the concentration mode can 
be much larger than that of the convective amplitude. 
The interaction is then dominated by the concentration 
mode. Here we describe this regime by considering the 
limit of weak diffusion of A, 77 <C 1. The fronts in A are 
then steep and the overlap between them can be ignored 
for any distance of 0(1). 

Based on numerical simulations |3l|j3^ ] and a prelimi- 
nary calculation for the interaction between widely sepa- 
rated fronts (77 = 0(1)) ^ we expect that stable pulses 
exist only if their velocity is opposite to the group veloc- 
ity s (for h > 0). Since the difference between these two 
velocities has to be small in a perturbative approach we 
are led to assume, in addition, that the group velocity s 
is small, s — rj 2 s 2 . Consequently, the resulting (average) 
pulse velocity v — (Otxr + 8txl)/2 will also be small. 
Here T — rj 2 t is a slow time scale and xr and xl are 
the position of the front at the right and at the left end 
of the pulse, respectively. Finally, in order to obtain a 
manageable expression for the solution of (@) we assume 



S = r? 4 (54, a = rj 2 a2, 



h = T] 3 h 3 , g = ri 2 g 2 . 



(9) 



This amounts to assuming the concentration mode to be 
small and its diffusion and decay to be slow as compared 
to the velocity of the front. These scalings lead to 

d t A + r] 2 s 2 d x A = tfdlA - A + cA 3 - A 5 + CA, (10) 
d t C = 7] 4 Sid 2 x C - V 2 a 2 C + r] 3 h 3 d x A 2 + rj 2 g 2 CA 2 . (11) 

Due to the smallness of the diffusion coefficient of A in- 
ternal layers arise at the transitions from the conductive 
to the convective state. Thus, in the presence of two 
fronts the system has to be separated into five regions. 
This is sketched in figj^. It shows also the typical shape 
of the fronts in the limit considered here. Whereas the 
traveling-wave amplitude is constant to lowest order in 
regions I, III and V, it changes rapidly from the con- 
ductive to the convective state in the transition regions 
II and IV. More precisely, in II and IV A and therefore 
also C varies on a fast space scale y = r\~ x x. In order 
to match the solutions in the different regions the width 
2A of regions II and IV needs to be small on the scale 

x but large on the scale y. This is achieved by taking 
A = ?/ 1 /2 Al/2 . 

The amplitudes A and C as well as the control param- 
eter c are now expanded as 



A = A + -qAi + 



c + rjci. 



(12) 



Here c is expanded around the value cq = 4/-\/3 for which 
in the absence of the concentration mode the conductive 
and the convective state coexist, i.e. neither state invades 
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the other. In the transition regions II and IV one obtains 
then at 0(1) 

= d*A + A o + c o A 3 -A%. (13) 

At 0(rj) and 0(r] 2 ) one finds, respectively, 

3 Tl A + s 2 d y A Q = (d 2 +£)Ai+ CiAo + Cl A 3 Q , (14) 



d Tl Ci = h 3 d y Al. 



(15) 



Here the intermediate time scale Ti = r/t and the lin- 
earized operator C = 1 + ScqAq — 5Aq have been intro- 
duced. Eq.((l|) yields A (y) = A f ((x R -x)/r]) in IV and 
Ao(y) = Af((x — xl)/t)) in II where the front solution is 

given by A f (() = A* (|(1 + tanhC)) 1/2 and A* 2 = y/Z. 
Its position is given by xr in IV and by Xl in II. Both 
can vary slowly in time, xr^ = xr,l(T). Due to trans- 
lational symmetry, the linearized operator d 2 + £ is sin- 



gular. Therefore eq.([14|) leads to a solvability condition 
in each of the two regions. Thus, by projecting eq.jl^) 
onto the left zero-eigenvector d y Ag one obtains evolution 
equations for xr and for xl, respectively. Note that in 
this approach the distance between the two fronts does 
not have to be assumed to be large in order to obtain 
two independent conditions. This is due to the small dif- 
fusion of A; any finite distance on the scale x is large on 
the relevant scale y. 

Since Aq depends on time only via xr^l (T) the general 
solution of eq. <RB) in II and IV is given by 



extremal values of the front solution At. Since C is non- 
zero in the regions in question, eq.(|18|) can always be 
solved for A -\ an d no solvability conditions arise in these 
regions. Eq.(|19|) has the general solution 



Ci = Ci(a;)e' 



-AT 



(21) 



with the amplitude C\(x) as yet undetermined. Note 
that in the present scaling the concentration mode de- 
cays only due to damping; the diffusive term is irrelevant 
despite the fact that C can vary on 0(1) length scales. 

To make use of the solvability conditions the integra- 
tion constants k^ 1 '^ appearing in the solution for C\ 
in II and IV are needed. They are obtained by match- 
ing the solutions in the different regions. In the limit of 
weak diffusion considered here < 1) the solution de- 
pends qualitatively on the direction of propagation of the 
two fronts. In the following we concentrate on the case 
in which the two fronts travel in the same direction and 
can therefore form a bound pair. If they travel towards 
each other no interesting interaction arises in this limit. 
For fronts traveling away from each other the calculation 
would have to be modified. For concreteness we assume 
both vr = OtXr and vl = QtXl to be positive. 

The traveling-wave amplitude vanishes in I and V and 
in III one has Aq = A*. Due to the damping term aC, 
the concentration mode vanishes also far ahead of the 
pulse, i.e. in V, . This implies k{ v = 0. At the right end 
of region III Ci is therefore given by 



Ci 



k{ 1JV (T) 



Otxr, 



A 2 



(16) 



with fc( /,/V (T) being integration constants in the respec- 
tive regions. Note that this solution is only valid as long 
as BtXr^l is O(l) or larger. For smaller velocities and in 
particular for stationary pulses the solution of the con- 
centration equation becomes more complicated, in par- 
ticular when the diffusion of C becomes relevant. 

In regions I, III and V the traveling-wave amplitude 
A is constant to lowest order. Therefore, the source 
term d x A 2 in the equation for C becomes less important 
than the damping term —aC. Thus, in these regions 
cqs. ([To|JTT| ) lead at 0(1), 0(r]) and 0{rf), respectively, 
to 



= A + c Al - Al, 
d Tl A = CAi + dA + ciAl, 
8tC\ — —aCi 



(17) 
(18) 
(19) 



where we have introduced the renormalizcd damping co- 
efficient 



a = a 2 - g 2 Al. 



(20) 



Note that for free-slip-permeable boundary conditions 
one obtains g 2 > p^ , i.e. the damping is reduced 
by convection. The solution of cq. (|l7|) is given by the 



C 1 {x R -A,T) = - 



vr(T) 



A* 



(22) 



To obtain C\ at the left end of region III at time t, i.e. 
at xl(T) + A, one has to take into account that the 
concentration field at that location has been generated 
by the leading edge of the pulse (which is now at xr(T)) 
at an earlier time T r = T — AT and has suffered an 
exponential decay since then. Thus one obtains 



Cx(x L + A,T) 



A* e~ &AT , (23) 



d T x R (T-AT) 
where the delay AT is determined by 



xr{T — AT) = xl(T). (24) 
Thus, in region II the concentration field C\ is given by 



Ci - h 3 A* 



1 



(1 



Al 



1 



iAT 



v L {TY A* 2 ' v R (T — AT) 



(25) 



For the solvability conditions C± need not be determined 
in region I. 

Projecting eq.(|lj) onto 8 v Aq leads in region II and IV, 
respectively, to the equations 
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v L (T) = s 2 
Vr(T) = s 2 



7 



27e" 



iAT 



v L (T) vr{T — AT) 



P, 



vr(T) 
d T L = v R (T) - v L (T) 

with L = xr — xl and 

7 = Vsh 3 , p = VSci. 

Here we have used 

(d y A ) 2 dy = A* 2 /4 = V3/4. 



(26) 

(27) 
(28) 

(29) 
(30) 



In terms of the front velocities VR^{t) the delay AT is 
given by 



/ v R {T')dT' =L{T). 
Jt-at 



(31) 



ThuSjin the general case, in which vr depends on time, 
eqs. (H© with © represent a set of implicit integral- 
delay equations. In the present scaling of weak diffusion, 
however, the dynamics in region IV is decoupled from 
that in region II. Therefore vr is constant and the re- 
tarded velocity vr(T— AT) is equal to vr. Consequently, 
the delay is given by AT = L/vr and one obtains 



vr = s 2 hp, 

v L [T) = s 2 — - H p, 



vl(T) 
d T L = v R - v L (T). 



vr 



(32) 

(33) 
(34) 



The results for negative front velocities are obtained in 
a similar way. Both can be combined into equations for 
the velocities v^t of the leading and the trailing front, 
respectively, 



vi = s 2 



7 

H 



+ sgn(vi)p 



(35) 



ry e -aL/\v t \ 

v t ■ = s 2 - — + 27 — |— j sgn(v t )p, (36) 



8tL 



1 



1 



7 27- 

Vt Vi' 



,-&L/\vi\ 



Vl 



2p. 



(37) 



Here we have used dxL = sgn(vi)(vi —v t ) which assumes 
that both fronts travel in the same direction. Note that 
eq.(|37|) does not represent a closed equation for L(T) 
due to the dependence of vt on L. The equation for v t 
is quadratic in vt and contains in addition the sign of v t , 
sgn(v t ). Elimination of v t in favor of L would therefore 
lead to complicated expressions involving the exponen- 
tial within square roots. We therefore prefer to discuss 
eqs.( |35|j3^ ) for the velocities instead of the equation for 
L(T). 



6||3_7j) 



It should be pointed out that the validity of eqs 
is not limited to large values of L. Therefore 7 e 
can be of the same order as 7 itself. In the present scal- 
ing the traveling-wave amplitudes A of the two fronts 
do not overlap on this length scale and do therefore not 
contribute to their interaction. We will comment on the 
possible consequences of an overlap in A and of dispersion 
at the end of sec.V. 



IV. DISCUSSION OF THE FRONT EQUATIONS 

In the limit L — > 00 eqs.( |35| , j36|) describe the dynamics 
of a single front. More precisely, a single front connecting 
the state A = A* for x — > — 00 with the state ^4 = 
for x — * +00 corresponds in the present calculation to a 
leading front for v > (cf. fig||a) and to a trailing front 
for v < (cf. figj|b). Its velocity v is therefore given by 
the relationship 



* 2 



or equivalently by 



(s 2 + p ± \J (s 2 + p) 2 - 4:.sgn{v)"fj 



(38) 



(39) 



Fig.[| gives a typical behavior of v as a function of p for 
such a front. It clearly exhibits a hysteretic transition 
from a 'fast' front to a 'slow' front when decreasing p 
from positive values through the saddle-node bifurcation 
at 



Psn = 2 



7 



s 2 . 



(40) 



It is expected that only fronts on the top and on the bot- 
tom branch are stable and those on the middle branch are 
unstable. Since the main difference between the fast and 
the slow front is the strength of the concentration mode 
the destabilizing mode of fronts on the middle branch 
will mainly involve the concentration mode. 

The possibility of a hysteretic transition between fast 
and slow fronts is not unexpected. The source of the con- 
centration field is given by the gradients of the traveling 
front in the traveling-wave amplitude. If the front is trav- 
eling fast, the source remains at any given location only 
briefly and therefore generates only a weak concentration 
field which in turn has only little effect on the propaga- 
tion velocity of the front. If the front is slow, however, 
a strong concentration field can be built up which slows 
down the front considerably providing a positive feed- 
back mechanism. Note that, in experiments, increasing 
p corresponds to increasing the Rayleigh number. Thus, 
in principle it may be possible to observe this drastic 
manifestation of the concentration mode experimentally. 

A right-traveling front connecting A — A* for x — ► — 00 
with A = for x — > +00 is not equivalent to one con- 
necting A — with A* traveling to the left (cf. fig.||c). 
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This is due to the fact that A is the amplitude of right- 
traveling convection rolls. Reflection in space therefore 
yields left-traveling rolls with opposite group velocity s 
and opposite advection of the concentration field, hd x A 2 . 
The velocity of such a front is obtained from eq.(|39|) by 
flipping the sign of p. The velocities for both types of 
fronts are shown in fig.|j. In contrast to fig.|] the velocity 
of the leading front is now marked by heavy solid lines 
and that of the trailing front - in the limit L — > oo - 
by thin dashed lines. For p = the convective and the 
conductive state have the same energy in this dispersion- 
less case. In the absence of any interaction the trailing 
and the leading front have therefore the same velocity. 
Due to reflection symmetry the concentration field af- 
fects the velocity of both types of fronts in the same way. 
Thus, vi{p) and vt(p, L — > oo) intersect at p = even in 
the presence of the concentration field, as indicated by 
the three open circles. These three solutions at infinite 
L imply the existence of three different localized waves 
(pulses) within this framework. 

At finite distances the fronts interact with each other. 
Fig.^| shows the somewhat complicated behavior arising 
from the interaction. The broken lines give the velocity 
Vt of the trailing front as a function of p and with the 
velocity vi determined correspondingly. The parameters 
are given by 7 = 0.15, a — 0.01, g — and s = 1. 
The distance between the fronts is L = 10 in figj^a and 
L = 30 in fig.^b. Eq.j36]) is quadratic in v t . One therefore 
obtains up to two solutions for given v\, p and L, i.e. the 
curve v t (p) can split up into two branches. Eqs.(j35 36) 
are, however, only valid if Vt and vi have the same sign. 
Thus, for v t < only one branch remains which is marked 
by heavy short dashes. Due to the interaction it is shifted 
to smaller values of p. Consequently, the intersection of vi 
and Vt, which signifies a pulse solution, occurs for p < 0. 
For vt > both branches are relevant. They are marked 
by heavy dashed and dotted-dashed lines, respectively. 
Due to the double-valued nature of vi one obtains four 
different values of Vt for a given value of p, which are 
associated in a pairwise manner with the upper and the 
lower branch of vi , respectively. Pulse solutions are given 
by the intersection points of the curves. Therefore, care 
has to be taken which value of vt is associated with which 
branch of vi. The relevant intersections are marked by 
solid circles. 

The dependence of the velocity curves on the length 
L allows conclusions about the character of the interac- 
tion. For pulse A in fig.^ it is repulsive; without inter- 
action (thin, dashed line) the magnitude of vt would be 
larger than that of Vi and the trailing front would catch 
up with the leading front. With decreasing distance L, 
however, the interaction strength increases and the mag- 
nitude of v t is reduced. Thus, the trailing front will ap- 
proach the leading front until their velocities match. If 
the distance were to be decreased further, the trailing 
front would start lagging behind until the equilibrium 
distance is reached again. This suggests that the repul- 
sive interaction leads to a stable pulse. By contrast, for 



pulse B the interaction is attractive and should render it 
unstable: the closer the trailing front comes to the lead- 
ing front the faster it travels. For pulse C the interaction 
is again repulsive. 

The above expectations about the stability of the 
pulses are borne out in a direct linear stability analysis. 
A linearization of eqs.(|35|-p7|) around the steady pulse 
solution with velocity vq and length Lq (and correspond- 
ing delay ATq = Lq/vq) yields for the growth rate a of 
perturbations 



2a 



1 - 



-&Lo/\vo\ 



(41) 



For a > the pulses are therefore stable as long as 
|^o It -1 < 10- The saddle-node bifurcation of the lead- 
ing pulse occurs at vo\vo\ — 7- This stability result cor- 
responds therefore to the two cases with repulsive in- 
teraction discussed above (A and C). Only the case 
vq < will lead to stable pulses since on the lower branch 
(0 < v \vq I < 7) the leading front itself is unstable to per- 
turbations in the concentration field which take it either 
to the fast or to the slow branch. This instability is not 
contained in the present analysis of the velocity and in- 
teraction of fronts (eqs.( |35| , |36| )), since it involves a change 
in the character of the front which is not merely given by 
a change in its position. 

Can one understand the stability of backward traveling 
pulses intuitively? Consider the experimentally relevant 
case 7 > . Then the negative gradient of A 2 at the 
right front leads to a negative peak of C whereas the left 
front generates a positive peak. If the pulse is stationary, 
both peaks have the same magnitude due to reflection 
symmetry as indicated in fig.|l| by the dotted line. At 
the right front the concentration field reduces the local 
growth rate of A and therefore pushes the pulse to the 
left. Similarly, the peak at the left front enhances the 
growth rate und also pushes the pulse to the left. Thus, in 
the stationary case both peaks change the pulse velocity 
by the same amount and therefore do not contribute to 
the interaction. If the pulse travels, however, this is not 
the case any more. For negative velocities both peaks 
become steeper to their left side and much shallower to 
the right (dashed line), and the concentration field of 
the left front influences the right front much more than 
the other way around. Thus, while the left front enters 
a region with C = 0, the right front propagates into a 
region with C > 0. The magnitude of the concentration 
field at the right front is therefore reduced as indicated 
by heavy lines in fig.0. Consequently, the right front 
is pushed less to the left than the left one resulting in 
a repulsive interaction. This argument shows that to 
lowest order the stability of the pulses requires that their 



2 Note that for a < the basic equations (|1| 
appropriate. 
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velocity v be opposite to the group velocity s: they travel 
backwards. 

The dependence of the length of the stable pulse on 
p is shown in figj|. As expected the length diverges at 
p = 0. It monotonically decreases with increasing \p\. 
Strikingly, within this framework the stable pulse exists 
for arbitrarily large negative values of p. To wit, for large 
p and small vi the velocity and length are given to leading 
order by 



vi = j/p + 



L = 



S2_l_ 

7a p 2 



(42) 



implying that the repulsive interaction becomes arbitrar- 
ily large in this limit. This is due to the prefactor 
in eq.(36): with decreasing magnitude of the velocity the 
build-up of the concentration field becomes stronger and 
stronger. By contrast, for the (attractive) interaction of 
fronts in the real Ginzburg-Landau equation alone the 
prefactor of the exponential is constant and the strength 
of the interaction remains bounded. 



V. NUMERICAL RESULTS 



In the derivation of (35 36|) we have pushed the diffu- 
sion of A as well as the diffusive term in C to higher order. 
To get an impression of the influence of these effects on 
single fronts and pulses we investigate eqs.(|^J|) numeri- 
cally. In fig.^ the results for the velocity v of a single front 
are given as a function of c = (4 + ryp)/\/3 for increasing 
values of 8 j4(| . The large symbols denote the velocity of 
a left-front, i.e. a front connecting A(x — > —00) = with 
A(x — > +00) = A* , whereas the small open circles denote 
that of a right-front. These results show that the singu- 
larity in c = c(v) at v = for a single front disappears 
when the diffusive term for C is included; instead the slow 
fronts on the bottom branch also undergo a saddle-node 
bifurcation and the two branches, which are separated in 
the approximation leading to eqs.( ^5| , ^6| ), are connected 
in an S-shaped curve. With increasing 8 the S-shape be- 
comes less pronounced and one obtains eventually a con- 
tinuous transition to a single-valued velocity. A similar 
transition occurs with decreasing coupling strength. 

The diffusive terms have also a strong effect on the 
regime of existence of the stable pulse. The results of 
numerical simulations of eq.((7|J|) are shown in fig.||(a,b) 
which give the length L and the velocity v of the pulse as 
a function of c for different values of 8. While the length 
appears to diverge with increasing c, the pulses disappear 
in a saddle-node bifurcation at small values of c. 

This behavior can be understood qualitatively in the 
following way. In the presence of diffusion not only the 
trailing but also the leading front is affected by the inter- 
action. For very weak diffusion, the effect on the leading 
front will, however, be negligible. Nevertheless, the solu- 
tion branches cannot be separated into trailing and lead- 
ing fronts any more. Instead one has to label them again 



as left and right fronts with velocities Vl and vr. Con- 
sidering the diffusionlcss result, we expect that for weak 
diffusion the stability of the pulses requires that their 
velocity be negative and that the isolated left-front be 
slower than the right-front. Therefore, c must be larger 
than a minimal value c m i n at which vl vanishes and be- 
low a maximal value c max at which vl = vr. This is 
sketched in fig.fi]. The difference Av in the velocities has 
to be compensated by the interaction term, which will 
now also depend on 8. Keeping the main ingredients of 
the interaction in mind one may expect it to be of the 
form f(v)e~ &L '' v ^ where f(v) remains now bounded for 
v — ► 0. Balancing the two terms one obtains the rela- 
tionship L = -cr 1 |i;|ln(Ai>//H). Thus, if Av -> at 
finite v the length diverges. This is the case at c max . 
If, however, v — > at finite Av, as is the case at c m i„, 
L -> 0. 

Of course, when L reaches 0(ri) the overlap in the 
traveling- wave amplitude A has to be taken into account, 
as well. Since for very small lengths all pulse-like initial 
conditions will collapse, a short, unstable pulse has to ex- 
ist which separates the basins of attraction of the pulse- 
less solution {A = 0) and that of the stable pulse. With 
decreasing c the stable pulse is therefore expected to col- 
lide with this unstable pulse in a saddle-node bifurcation 
for c slightly above c m i n ■ The unstable pulse should cor- 
respond to that of the uncoupled real Ginzburg-Landau 
equation where it exists for c > c ma x and has diverging 
length for c — > c+ . In the presence of the concentration 
mode it appears to be created together with the stable 
pulse already for c < c max at a finite length. 

With increasing 8 - and also with decreasing h - the 
hysteresis in the front velocities decreases. The range 
(cmimCmax) will therefore decrease and vanish above 
some critical value of 8. Thus, stable pulses exist only 
for sufficiently strong coupling h and sufficiently weak 
diffusion 8 (cf. fig.||)). 

In principle, it is conceivable that the left-front solution 
reaches its saddle-node bifurcation at negative v^. In 
that case the branch of stable pulses has to cease because 
no suitable leading front exists. The length at which 
this occurs could be sufficiently large that the overlap 
in A is irrelevant and may diverge with increasing 8 or 
decreasing h. 



VI. CONCLUSION 

In this paper we have investigated the dynamics of 
fronts which interact via a long-wavelength mode. The 
work was motivated by convection in binary liquid mix- 
tures where the fronts connect the motionless state with 
the convective state and the long-wavelength mode corre- 
sponds to a concentration mode which arises due to small 
mass diffusion. We have shown that this additional mode 
can lead to stable pulses in the form of bound pairs of 
fronts even in the absence of dispersion. 
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How do these results relate to experiments in binary- 
mixture convection? Of course, the physical system is not 
dispersionless, i.e. the coefficients in equation ([!]) as well 
as the traveling-wave amplitude are in general complex. 
It has been shown previously p9pc| ] that this by itself can 
lead to stable bound pairs of fronts without the presence 
of a concentration mode. For weak dispersion one obtains 
an equation of the form [29|j3(| 

d T £ = fci(c-c )-fc2e- L/ « + ^ (43) 

for the length L of a pulse with k\ and ki positive. The 
exponential term represents the attractive interaction be- 
tween fronts due to the overlap of their exponential tails 
in the convective amplitude, while the algebraic interac- 
tion term arises from dispersion. For suitable parameters 
it can be repulsive (£3 > 0). Within this framework, for 
/C3 > and c < cq one has two pulse solutions as in the 
dispersionless case discussed above. Again the one with 
larger L is the stable oneQ. This ordering appears to per- 
sist for strong dispersion where the stable pulse may be 
identified with a dissipatively perturbed soliton solution 
of the nonlinear Schrodinger equation pi] . 

In recent experiments, both stable as well as unstable 
pulses have been investigated in great detail pi] . The 
unstable ones have been obtained by employing a suitable 
servo control. In that work two regimes have been found. 
For a value of the separation ratio of ip = —0.127 two 
pulses exist for a given Rayleigh number and the stable 
pulse is the shorter one. For ip < —0.17 the branch with 
the longer, unstable pulses turns around and one is left 
with a single pulse solution which is stable. 

It is tempting to speculate that the experimental find- 
ing of longer, unstable pulses can be described by in- 
cluding dispersion as well as the concentration mode. 
The derivation of the relevant equations appears to be 
involved. Considering eqs.(|43"|) and ( |33 , [mB7]) one may 
expect that an equation of the form 

d T L = kx{c- c ) - k 2 e~ L ^ + — - 2le~ &L/M (44) 

L vi 

could capture the main aspects. With 7 > as before 
this equation allows four solutions of forward traveling 
pulses (vi > 0). This is illustrated in fig|| where the 
dL/dt is plotted as a function of L for two exemplary 
cases. Again the shortest pulse at L = 0.25 (marked by 
an open triangle) is unphysical. The unstable pulse A 



constitutes the separatrix between the pulseless conduc- 
tive state and the state with a stable pulse (B). It has 
not been observed. We identify the experimentally ob- 
served pulses with the stable pulse B and the unstable 
pulse with pulse C. Thus, within the extended equations 
an unstable longer pulse can coexist with a shorter stable 
one as seen in the experiment. Reducing c slightly below 
Co leads to a second stable pulse (D). Since its length 
diverges as c — > Co it may be too long to be observed in 
the present experiments. 

When j/vi > is decreased towards the two pulses 
B and C merge in a saddle-node bifurcation and - in 
particular for backward traveling pulses (v < 0) - one is 
left with the unstable pulse A and the stable pulse D. 
We identify this situation with that found in the exper- 
iments for ip < —0.17. This interpretation is supported 
by the fact that in the experiments the stability of the 
long pulses changes very close to where their velocity 
changes sign. Within the present scenario one would ex- 
pect the length scale for the long, stable pulses to be 
set by the decay length £ = \v\/dt2 of the concentration 
field. Since convection has a strong influence on the decay 
rate through a = ct2 ~ 92\Ao\ 2 , such an estimate requires 
more than the knowledge of the Lewis number and the 
measured velocity of the pulses. It is therefore not too 
surprising that taking a-i alone instead of a leads to an 
unreasonably short decay length I. Clearly more detailed 
calculations are necessary to substantiate the discussion 
based on eq.@ @. 

In this paper we have concentrated on localized solu- 
tions which can be considered as bound pairs of fronts, 
each of which forms a heteroclinic orbit in space. This 
requires, of course, a subcritical bifurcation. It is 
worth noting that numerical simulations of the extended 
Ginzburg-Landau equations (^,^|) have shown that local- 
ized waves can also arise in the case of a supercritical 
bifurcation, for which no single fronts exist J32J. These 
pulses correspond to homoclinic orbits in space. So far, 
they have not been captured in any analytic description. 
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3 In the regime in question eq.(^3|) has in fact three solutions. 
The third solution, which is formally stable within eq. (fig]) , is, 
however, always confined to L/£ = O(l). This small distance 
between the fronts is inconsistent with the assumption of weak 
interaction between them since &2 = 0(1). In fact, this so- 
lution is unphysical; as discussed in the previous section, the 
smallest pulse must be unstable. 
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FIG. 1. Influence of the propagation of the pulse on the 
concentration field. A stationary pulse has symmetric con- 
centration field (dotted line for C). In a propagating pulse 
(dashed line for C) the trailing peak of C is reduced leading 
to an effective repulsion between the fronts (see text). 

FIG. 2. Sketches of leading and trailing fronts for v > 
and v < 0. a) leading front with v > 0. b) trailing front with 
v < 0. c) trailing front with v > 0. d) leading front with 
v < 0. 

ffi 



FIG. 3. Dependence of the velocity of a single front on p. 

FIG. 4. Velocity of a leading (solid line) and of a trailing 
front (dashed lines) as a function of p. Thin dashed lines 
denote vt for L — » oo. Thick dashed lines give v t for L = 10 
and L = 30, respectively. 



FIG. 5. Dependence of the length of pulses on p. 

FIG. 6. Numerical results for the velocity v of a single front 
on c. The parameter values are s = 0.3, rf — 0.05, a — 0.01 
and h = 0.03. 

FIG. 7. Sketch of the influence of diffusion and damping 
of C on the velocities vi, t of fronts. Solid lines for L — > oo, 
dashed line gives vt for finite L. 

FIG. 8. Numerical simulations of pulses with s = 0.3, 
if = 0.05, a = 0.01 and h = 0.03. The grid size is Ax = 0.1 
and the time step is At = 1. The simulations ar done in a 
frame moving with the pulse. 

FIG. 9. dL/dt as a function of L according to eq.(Q). In 
this regime eq.(E3) has up to 4 steady pulse solutions. 
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